Conservation Genetics

EU has reclassified the sika deer (Cervus nippon) as an undesirable invasive species based on reports that hybridization with the indigenous red deer (C. elaphus) may produce fertile offspring. Since sika-derived DNA previosuly introduced into the red deer population (introgression) cannot be removed, the crucial question is whether new (F1) hybridisation occur. To address this, we analysed the chromosomes in 56 sika and 22 red deer. All red deer had a chromosome number 2n=68. In contrast, the chromosome number in sika ranged from 64 to 67, due to the variable presence of two sika-specific Robertsonian translocations (ROB1,ROB2). In the free-ranging sika population in Jutland, >90% of the sika deer were homozygote for at least one of these ROBs, excluding that they could be F1-hybrids. Moreover, ROB2 was in Hardy-Weinberg equilibrium, further supporting the absence of gene flow between the two species. In contrast, ROB1 was in Hardy-Weinberg disequilibrium, suggesting negative fitness of heterozygotes, including potential F1-hybrids. In Jaegersborg Deer Park, the eight examined sika deer had the same genotype (absence of ROB1, homozygosity of ROB2), supporting that it is a founder population which may have been isolated for [~]100 years. Again, none of these can be F1-hybrids due to the homozygosity of ROB2. We conclude that F1-hybridisation between sika and red deer either does not occur or occur very rarely in Denmark. The study establish the Danish sika-populations as unique models for adressing important biological questions: What underlies the absence of hybridisation? Why are ROBs frequent in sika deer but not in the closely related red deer? How fast do new species/subspecies develop in isolated founder populations? Which factors determine, that some ROBs have little heterozygous effects, whereas others are selected against, with implications for the role of ROBs as genetic barriers promoting speciation, and for fertility problems in some human ROB carriers.

Rosa gallica L., the French rose, is a perennial, tetraploid, heterozygous species that naturally propagates by seed and sucker. It occurs in the wild, primarily in Europe, and also exists as cultivated varieties. R. gallica cultivars were extensively bred and cultivated in France at the beginning of the 19th century. Although several hypotheses have been proposed regarding the species expansion based on historical records, none have been assessed using molecular data. Indeed, its genetic diversity has so far been investigated only at local or regional scales, hindering the identification of the evolutionary factors shaping its present-day distribution. Using 29 sequenced microsatellites, we genotyped a comprehensive sample of 1618 individuals, including wild R. gallica from 219 sites across the species range, rose cultivars, and specimens from other Rosa species. We then detected clonal lineages and characterized the range-wide genetic diversity and structure, aiming to disentangle the roles of natural and human factors in shaping the distribution of R. gallica, with particular focus on France. French diversity appears particularly structured compared to the rest of the range, suggesting multiple origins within France. Populations in South Alps, Central Eastern Europe, and Eastern France appear to have recolonized naturally from a single southern glacial refugium. In contrast, populations in the western part of France likely resulted from more recent natural or human-mediated dispersal. Finally, clonal lineages containing both wild and cultivated individuals were predominantly found in France, highlighting the role of human-mediated dispersal in 28 of the 98 French sites studied. These findings show that the present-day natural range of R. gallica was shaped primarily by post-glacial recolonization, but also reveal a contribution of human activities to its recent dispersal, particularly in France, where cultivated varieties were intensively bred and exchanged.

Understanding how contemporary population declines affect the genomic diversity and structure of threatened species is important for effective conservation. The Eurasian curlew (Numenius arquata) is experiencing severe population declines across Europe, with Ireland among the most extreme, showing declines exceeding 90% over 40 years. Genomic data are increasingly incorporated into policy and used to assess conservation status by estimating genetic diversity, differentiation, inbreeding, effective population size, and adaptive divergence. Such data for curlew is scarce, and the population structure among northern and north-western European breeding populations remains unclear. To address this, we generated whole-genome resequencing data for 56 curlews across Ireland, Britain and Sweden. Irish and British populations showed minimal interpopulation differentiation, but both were substantially differentiated from Sweden. This was apparent from principal component analysis, and admixture and FST analyses. Measures of genetic diversity (nucleotide diversity, heterozygosity, Watterson's{theta} ) were similar across populations. A slightly elevated Tajima's D in Ireland, along with elevated FROH in Ireland and Britain relative to Sweden, may be the early genomic signs of recent population declines. We identified locally selected candidate genes. These had putative roles in metabolic processes, the immune response, and were potentially associated with distinct migratory behaviours and environmental conditions. We find a potential lag in genomic effects of decline being detectable following population contraction. We also show highly migratory species can exhibit differentiation in ecologically relevant traits, potentially driven by high site fidelity. These findings warrant consideration in translocation planning and broader conservation strategies.

Due to the sparse and uneven availability of genomic resources, it remains challenging to appraise genomic attributes for species of conservation concern. While long-read sequencing enables assessment of genetic diversity at unprecedented scale, accessing high-quality tissues remains a challenge for non-model species. Here, we explore an alternative sampling strategy for tissues where "gold-standard" cryopreservation is infeasible. Focusing on the Australian scincid lizard Cryptoblepharus pulcher, we compare DNA obtained from various ethanol-preserved tissue types and DNA extraction kits, and ask whether high-molecular weight DNA can still be retrieved. PacBio HiFi sequencing of the most promising sample yielded a highly contiguous reference-level assembly, validating this approach in vertebrates, specifically lizards. After scaffolding the assembly to chromosome-level, we then used a comparative approach to shed light on the evolution and demise of C. egeriae, a closely related, now Extinct-in-the-Wild species. Surprisingly, despite being a wide-spread continental analogue with a similar ecology, C. pulcher has lower genetic diversity and long-term historical population size than C. egeriae, an island endemic. However, C. pulcher also shows fewer runs-of-homozygosity, supporting prior reports that C. egeriae experienced recent inbreeding. Together, these findings demonstrate that a practical and cost-effective preservation strategy can still yield high-quality genomic resources in vertebrates, as well as valuable insights that are relevant in an age of biodiversity decline.

Among Madagascar primates, the sportive lemurs (family Lepilemuridae) have seen their species diversity increase from eight in 2005 to 26 in 2009 mostly by applying the phylogenetic species concept to DNA barcode data. Despite the genus being speciose, only one case of sympatry is known from northern Madagascar, where two sportive lemur species described based on low mtDNA divergence, Lepilemur ankaranensis and Lepilemur milanoii, were found to co-occur at the center of their joint distribution range. Here, to clarify the taxonomy of these two species and examine their sympatry, we apply an integrative taxonomic framework to genomic and morphological data from 84 individuals of L. ankaranensis and L. milanoii, encompassing their entire distribution range and the forest of Analafiana, beyond their southernmost limit. Using clustering, multivariate, and isolation by distance analyses, we find no evidence of a sympatric zone and show that despite clear genetic differentiation between regions, the genomic and morphological diversity of the L. ankaranensis, L. milanoii-Analafiana group is clinal and explained by geographic distance. These results clarify that L. milanoii is a junior synonym of L. ankaranensis and that the Analafiana forest population belongs to L. ankaranensis, extending its distribution. It further implies that the 'sympatric' zone, the Andrafiamena forest, hosts conspecific individuals with slightly differentiated mtDNA backgrounds, rather than slightly differentiated sympatric species. Lastly, we re-evaluate the IUCN conservation metrics of L. ankaranensis, which continue to qualify as Endangered (EN) under the B1ab(i-v) criteria.

High-throughput amplicon sequencing has transformed microsatellite (STR) genotyping by overcoming many of the limitations of fragment-length analysis, enabling more accurate, cost-effective, and standardized genotyping. Yet, protocols specifically designed for high-throughput sequencing (HTS)-based STR genotyping from low-template and degraded DNA remain scarce, despite the prevalence of these challenging sample types in ecological and conservation contexts. We present a methodology for the de novo development of robust STR multiplex panels together with a laboratory protocol for efficient and reliable STR genotyping by sequencing with low quantity and quality DNA samples. The protocol comprises (i) an automated bioinformatic pipeline to design large sets of short tetranucleotide markers optimized for multiplex amplicon sequencing of degraded and low-template DNA; (ii) guidelines for efficient in vitro optimization of multiplex amplification using directly low quantity/quality template DNA; and (iii) a library preparation procedure that improves detection of low-level allele signal while enabling quality assessment of STR amplicon sequencing under limiting DNA conditions. We demonstrate the approach by developing and validating STR panels for non-invasive genotyping of three large carnivore species: a 44-plex for the grey wolf (Canis lupus), a 41-plex for the Eurasian lynx (Lynx lynx), and a 30-plex for the brown bear (Ursus arctos). Multiplex performance was high, with [≥]91% of samples successfully genotyped at [≥]50% of loci (allele size range 28-110 bp across panels) and correctly assigned to known individuals, negligible levels of noise in the controls, and high discriminatory power (PIDsibs [≤]2.4 x 1e-12), also owing to sequence variation among same-length alleles at 15-50% of loci. The approach is broadly applicable to animal and plant species, a wide range of sample types, and large-scale analysis such as genetic monitoring. Our study reinforces the value of STR amplicon sequencing for ecological and conservation applications while highlighting the importance of marker design and laboratory workflows tailored to HTS-based genotyping for accurate and efficient implementation.

Genomic signatures can provide key insight into the evolutionary history and remaining adaptive potential of threatened populations. As demographic decline erodes both diversity and the processes maintaining it, understanding how remaining variation is distributed becomes increasingly important for conserving species like the staghorn coral, Acropora cervicornis, a foundational but critically endangered Caribbean reef-builder. We analyzed 46 high-coverage A. cervicornis genomes from 10 locations across the tropical western Atlantic to evaluate neutral and adaptive structure, genomic diversity, demographic history, inbreeding, and connectivity, and generated a regional haplotype reference panel for future genomic monitoring. Genome-wide analyses recovered recurring regional substructure, but differentiation was modest and partly explained by isolation-by-distance and spatial variation in effective migration. Subpopulations had similar levels of genomic diversity, shared demographic history, and limited evidence of local adaptation. These patterns support interpreting sampled Caribbean populations as a single evolutionarily significant unit (ESU) containing multiple regional management units (MUs), rather than as deeply divergent evolutionary lineages. Despite substantial retained variation and low current inbreeding, estimated contemporary effective population size was small, suggesting an increased vulnerability to the effects of drift as demographic collapse continues, especially if structure is reinforced by isolated management. Together, our findings emphasize the urgent need for interventions that preserve and enhance genomic diversity, including risk-managed assisted gene flow. Supported by the haplotype reference panel developed here, these strategies will require coordinated efforts across regional entities to conserve and restore A. cervicornis as a jointly managed, single ESU.

The red fox, Vulpes vulpes, is abundant in England and can exert limiting effects on their avian and mammalian prey. Large-scale gamebird releases in the UK may be sustaining high predator numbers, leading to greater predation on other prey species, especially once gamebird stocks are depleted. However, little is known how gamebird release affects the broader diet of the red fox. Here we used DNA metabarcoding to assess the diet of foxes from 18 agricultural estates in Southern England, 10 of which released large numbers of gamebirds (red-legged partridges - Alectoris rufus, and common pheasants - Phasianus colchicus) and 8 which did not. Scats were collected over one year, allowing for seasonal investigation of the foxes' diet. We investigated the vertebrate species consumed and compared the non-gamebird dietary diversity and composition between release and non-release estates and across seasons. The field vole (Microtus agrestis) was found to be the most frequently predated species overall. Brown hares and field voles were detected more on release sites, while bank voles and dog faeces were detected more on non-release sites. We found little evidence that foxes predate ground nesting birds or other species of concern. The dietary diversity was significantly lower on estates that released gamebirds, and this difference was most notable during the post-shoot, spring months (February to April). On both estate types, diversity was highest in the summer months. The altered predatory behaviour due to gamebird release is likely to affect the populations of the non-game prey of the red fox and this should be considered when designing policies regarding gamebird management and biodiversity conservation. Alternative predation combined with predator control may be reducing predation pressure on non-game prey where gamebirds are released, but high fox density elsewhere likely results in higher predation pressure on a wide range of species.

Understanding the evolutionary history of species within a geographic context is key to historical biogeography, as it reveals how geological and climatic changes shaped biodiversity. This is especially important in ecologically significant regions like the Amazon and Orinoco basins. Together, they host the worlds greatest freshwater fish diversity ([~]3,500 species), sharing a common but not yet fully understood evolutionary history. The gilded catfish (Brachyplatystoma rousseauxii), an ancient species widely distributed as a metapopulation in Neotropics, is an important model for studying past connectivity, divergence, and historical processes shaping fish diversity between these basins. This study analyzed the genetic structure, connectivity routes, and demographic history of B. rousseauxii using nuclear (microsatellite and ddRADseq) and mitochondrial DNA. Population structure analyses and coalescent models indicate that B. rousseauxii populations from the Orinoco and Amazon basins are genetically distinct, with no evidence of current gene flow. However, our results support the occurrence of a possible secondary contact event after the divergence, with the Boa Vista population retaining the genetic signal of this process. The ancestral population split occurred at the Rupununi Portal around 2.54 Ma (ddRAD) or 1.31 Ma (mtDNA). Then, the species colonized the Branco and Orinoco Rivers [~]1.90 Ma (ddRAD) or 0.6 Ma (mtDNA), rapidly expanding in the Orinoco (>1.3 or >0.29 Ma), while colonization of the Amazon from the Branco River was more recent ([≤]1.0 or [≤]0.15 Ma). Population expansion signal was detected in the Orinoco ([~]0.20 Ma), whereas the Amazon remained stable. Our findings suggest that the rise of the Vaupes Arch in the Late Miocene does not explain the observed genetic divergence. Likewise, the Casiquiare Canal and Japura-Guaviare headwaters are not connectivity routes between basins. Instead, the Rupununi Portal, including the recent capture of the Branco River by the Negro River, was the last point of connection and played a key role in shaping B. rousseauxiis distribution. These findings provide insights into Neotropical fish biogeography and the historical configuration of the Orinoco and Amazon basins.

Understanding the spatial patterns of biodiversity and their environmental determinants is fundamental to ecology and conservation, particularly in mountainous transitional zones where species assemblages can shift rapidly. Amphibians, as highly sensitive organisms, are excellent indicators of environmental change, yet quantitative assessments of their diversity in the Huangshan Mountains-a biodiversity hotspot at the junction of the Palaearctic and Oriental realms in southern Anhui, China-remain limited. We conducted systematic line-transect surveys across 200 grids (5 km x 5 km) during the spring and autumn of 2023 and 2024, recording a total of 10,342 individuals belonging to 23 species, 20 genera, 9 families, and 2 orders. Three species (Fejervarya multistriata, Microhyla fissipes, and Bufo gargarizans) were identified as dominant, and two nationally protected species (Hoplobatrachus chinensis and Andrias davidianus) were detected. Inter-annual alpha diversity did not differ significantly, but pronounced seasonal variation was observed, with spring supporting higher Shannon-Wiener and Pielou evenness than autumn in both years. Using generalized additive models and model averaging, we found that annual precipitation and the normalized difference vegetation index (NDVI) were consistently the strongest positive predictors of Shannon-Wiener diversity, Simpson dominance, and species richness, while species richness also declined significantly with increasing human footprint. Total beta diversity (Soerensen dissimilarity) was very high (0.981) and overwhelmingly driven by species turnover (97.86%) rather than nestedness. Partial Mantel tests and distance-based redundancy analysis further revealed that environmental distances-particularly annual precipitation-significantly shaped overall beta diversity and its turnover component after accounting for geographic distance. These results highlight the predominant role of climatic and vegetation gradients in structuring amphibian assemblages in this subtropical mountainous region, providing baseline data to inform local conservation strategies.

Heliococcus summervillei is an emerging invasive mealybug that causes severe dieback in grasses in pastures and turfgrass landscapes. It is widespread in Australia and has recently been detected across the Caribbean, Mexico, and the United States. Accurate identification of mealybugs is challenging due to cryptic morphology, overlapping diagnostic characters, and limited taxonomic expertise and literature, which makes molecular tools essential for regulatory diagnostics and management. We developed the first Cytochrome Oxidase I (COI) barcode for H. summervillei and used it to examine mitochondrial variation across available populations. COI sequences reveal approximately a 10.2% mitochondrial split between the Type A and Type B variants. Phylogenetic, haplotype network, and genetic distance analyses show that all invasive range populations share one haplotype associated with a recent invasion in the United States, Australia, Pakistan, and the Caribbean, whereas the Barbados lineage contains two closely related haplotypes that represent a historically stable mitochondrial variant. Together, these results establish the first COI reference library for H. summervillei, clarify mitochondrial lineage structure, and provide a practical barcode tool that enables rapid identification of invasive populations and supports timely regulatory and pest management responses. Recognizing mitochondrial variants also establishes a framework for resolving lineage-specific biological and management traits and strengthens reconstruction of introduction pathways central to regulatory decision-making and limiting further spread.

Predation is a major selective force shaping lemur behaviour and population dynamics, yet direct observations remain rare, particularly for small nocturnal species. Here, we document predation on mouse lemurs (Microcebus spp., Cheirogaleidae) in Andohahela National Park, south-eastern Madagascar, and using complementary evidence from radio-telemetry and owl pellet analyses. Two radio-collared individuals were confirmed as prey of endemic snakes, Ithycyphus oursi and Madagascarophis meridionalis (both Pseudoxyrhophiidae), representing the first documented records of Microcebus predation by these species. Examination of owl pellets and prey remains from three sites and from three owl species revealed a single predation event by Asio madagascariensis (Strigidae), likely involving M. tanosi, no evidence of mouse lemur predation by Tyto alba (Tytonidae) despite high local prey availability, and a single and first predation record from Athene superciliaris (Strigidae). Material from A. superciliaris roosts was otherwise dominated by invertebrates, indicating that primate predation is likely opportunistic. Together, these findings expand the known predator guild of mouse lemurs and suggest that snake predation is rarely detected, and its contribution to mouse lemur mortality and population dynamics is likely underestimated. Our findings demonstrate how combining behavioural field observations with dietary evidence can uncover otherwise undetected predation events and clarify predator prey relationships in nocturnal primates.

The White-bellied Heron (Ardea insignis) is one of the worlds rarest birds, with fewer than 60 known individuals remaining in the wild. Whether this extreme rarity reflects a recent anthropogenic collapse or a long history of persistently small population size has remained unknown, limiting our understanding of the species evolutionary resilience and conservation needs. Here, we present the first high-quality reference genome for A. insignis, generated using Oxford Nanopore long-read sequencing and complemented with Illumina whole-genome data. Comparative mitochondrial and nuclear phylogenomic analyses consistently recover A. insignis as the sister species of Purple Heron (A. purpurea), while revealing moderate mitonuclear discordance among deeper ardeid lineages. Genome-wide analyses demonstrate exceptionally low heterozygosity and extensive runs of homozygosity relative to the widespread and closely related Great Blue Heron (A. herodias), indicating pronounced genomic erosion and long-term inbreeding. However, the predominance of short and intermediate-length homozygous tracts, together with robust Pairwise Sequentially Markovian Coalescent (PSMC) reconstructions across alternative parameterizations, indicates that A. insignis has persisted with comparatively small effective population sizes over much of its evolutionary history rather than experiencing only a recent demographic collapse. The two sampled individuals nevertheless differ in the abundance of longer homozygous tracts, indicating that inbreeding accumulated over the past few generations has not been uniform among the surviving birds, despite their shared history of chronic rarity. Our results indicate that the White-bellied Heron represents a lineage that has survived prolonged demographic adversity and that its greatest genetic challenge may be limited adaptive potential rather than recent genomic deterioration alone. Beyond providing the first genomic resource for this critically endangered species, our study establishes an evolutionary baseline for future monitoring and highlights the importance of integrating genomic and ecological data to guide conservation strategies for species persisting at the edge of extinction.

Context Pumping stations pose a threat to fish globally through land use change, habitat fragmentation and entrainment risk, with the catadromous and critically endangered European eel particularly impacted. Objectives/methods Establish, model, assess and understand the present-day distribution of European eel and resident fishes in 152 pumping station catchments in a once extensive wetland (The Fens) using eDNA metabarcoding (855 samples over two and half years), with specific focus on anthropogenic influences on hydrological connectivity and habitat quality. A removal survey design maximised confidence in negative results while minimising time and consumable costs. Results Eel occurrence upstream of pumping stations was low (occupancy = 28.3%) and positively associated with catchment area, fish species richness and natural hydrological connectivity (gravity drainage or flooding) and negatively associated with distance from the tidal limit. Fish species richness replaced catchment area and improved model performance, potentially acting as a biotic indicator of habitat quality and connectivity. Pumped catchments with manually operated upstream water transfers had reduced eel presence, potentially linked to the direction of water flow or the timing of operation. By contrast, fish species richness increased in these catchments during summer, suggesting displacement into unsuitable long-term habitats. Physical habitat maintenance had no detectable effect on eel occurrence or fish species richness. Conclusions This study provides the first landscape-scale assessment of European eel distribution and drivers of occurrence in pumped river catchments. The highly novel and comprehensive insights have implications for European eel conservation as well as infrastructure and catchment management, including compliance with legislation (EC Regulation No. 1100/2007).

Population bottlenecks reduce genetic diversity, increase the fixation of deleterious mutations and elevate extinction risk. Identifying lineages experiencing bottlenecks is a key goal of conservation genetics, facilitating the allocation of limited resources to at-risk species. Although whole-genome sequencing has improved bottleneck detection by reconstructing demographic history, these methods often require extensive population sampling, limiting their application. Previous studies of species showing population bottlenecks have reported an increased number of pseudogenes in the olfactory receptor (OR) gene family, however whether such evolutionary dynamics can be used as comparative biomarkers of genomic decline remains unknown. By quantifying the number of lineage-specific duplication and pseudogenization events, we introduce the duplication-to-loss ratio (DLR), a comparative metric exploring the rate at which chemosensory gene loss is offset by the generation of novel receptors. We characterize the chemosensory repertoires of 21 felid species, including species with known historical bottlenecks, to establish the utility of this DLR metric. Subsequently, we evaluate its usage across additional mammalian families, specifically Ursidae and Pinnipedia, to determine its utility beyond Felidae. Our DLR metric recovers several felid species with a history of genomic decline, including cheetah (Acinonyx jubatus) and black-footed cat (Felis nigripes), as well as the giant panda (Ailuropoda melanoleuca), polar bear (Ursus maritimus), Hawaiian monk seal (Neomonachus schauinslandi) and northern elephant seal (Mirounga angustirostris). Our results demonstrate the utility of the OR gene repertoire as a scalable, robust biomarker for identifying comparative population decline, prioritising species for conservation genomic investigation using only the reference genome.

Understanding the spatial distribution of species and their patterns of endemism is necessary for establishing effective conservation priorities. Despite the vital pollination services bees provide, Californias native bee distribution patterns remain largely unexplored relative to plants and butterflies. We analyzed bee species richness and endemism across California and their concordance with plant distributions. Richness was high across areas of the California Floristic Province, including the Sierra Nevada, San Francisco Bay Area and Central Coast, South Coast Ranges, and the Transverse and Peninsular ranges. Bee endemism was more localized, concentrated in the San Joaquin Valley, eastern Sierra Nevada and adjacent Great Basin, Sierra Nevada foothills, and California deserts. Because richness and endemism appear to operate at different spatial scales and likely respond to different environmental drivers, effective conservation strategies must address both. Additionally, conservation plans that incorporate both plant and bee diversity are needed to achieve more comprehensive biodiversity protection.

Giant clams are some of the largest bivalve molluscs. They form a vital partnership with Symbiodiniaceae dinoflagellates that supply most of their energetic requirements. However, the factors that shape giant clam-associated photosymbiont communities remain unknown. Here, we profiled Symbiodiniaceae communities using ITS2 metabarcoding in eight giant clam species (Hippopus hippopus, H. porcellanus, Tridacna crocea, T. derasa, T. gigas, T. maxima, T. noae and T. squamosa) from 11 sites across the Philippine archipelago. Symbiodiniaceae community structure was shaped by an interplay between giant clam host and environment. Most giant clams were dominated by members of a single symbiont genus, with Cladocopium as the most prevalent, followed by Durusdinium and Symbiodinium. However, giant clam hosts also exhibited flexibility in their symbiotic partners that was evident across sites. Differences in giant clam-associated symbiont communities may contribute to differences in holobiont function and adaptability to variable environments. These findings deepen our understanding of giant clam-Symbiodiniaceae associations, offering a framework for predicting how giant clams may be affected by increasingly stressful reef conditions and, more importantly, informing strategies to improve mariculture and conservation practices.

The accelerated loss of biodiversity in sub-Saharan Africa threatens the functioning of tropical ecosystems. In Cote dIvoire, the Comoe National Park (PNCOMOE), a Sudano-Guinean savannah, and the Tai National Park (PNTAI), a dense rainforest, both UNESCO World Heritage Sites, are home to fauna assemblages of global importance, whose long-term resilience remains insufficiently quantified. This study assesses and compares, over a decade (2014-2025), the functional stability of vertebrate communities in these two contrasting ecosystems, using nine metrics covering resistance, invariance, persistence, interspecific synchrony, Tilmans stability, Jacobian resilience and a Composite Stability Index (CSI). Abundance data for 107 vertebrate species were collected via foot transects at PNTAI and aerial surveys at PNCOMOE. The stability metrics were calculated using the R package estar, integrated with an alpha diversity analysis (Shannon H', species richness S, Pielous evenness J') and a Jacobian spectral analysis within a multidimensional ecological assessment. PNTAI (0.708) exhibits significantly higher alpha diversity (H' = 2.82; S = 55.7 taxa) and community resilience 4.6 times higher than in the PNCOMOE (0.153). Its interspecific asynchrony index (0.504) reveals a strong portfolio effect, absent in the PNCOMOE (0.232). In contrast, PNCOMOE exhibits higher temporal invariance (0.382 versus 0.116) and Tilman stability (0.276 versus 0.152), reflecting more predictable dynamics. The overall ICS favours the PNTAI (0.484) and (0.370). The Jacobian analysis detects local instability in both parks (Re({lambda}max) = 5.58 at the PNTAI; 3.73 at the PNCOMOE). The two parks exhibit distinct yet complementary stability architectures: PNTAI relies on dynamic stability based on resilience and interspecific compensation, whilst PNCOMOE demonstrates conservative stability through temporal regularity. The absence of calculable resilience at PNCOMOE suggests a potential crossing of a functional degradation threshold, arguing for urgent restoration interventions and differentiated conservation strategies, tailored to the resilience mechanisms specific to each ecosystem.

Reliable population estimates are essential for wildlife management, yet monitoring schemes often do not match the administrative scale at which decisions are made. We illustrate this challenge using the Eurasian lynx in the canton of Valais, Switzerland, where official state monitoring is fragmented across three reference areas surveyed in different years. We analyzed three winters (2023-2026) of independent, canton-wide camera trap data, recording 899 independent lynx captures (27, 31 and 34 adults per winter). Lynx distribution and reproduction concentrated in the Northwest and connected to the thriving Pre-alpine populations. Density modelling for 2025/2026 estimated 37 independent lynx (95% CI: 26-52) on the whole cantonal territory, corresponding to a density of 1.09 (0.77-1.56) individuals per 100 km2. These estimates are substantially below the figures improperly extrapolated from a cross-cantonal reference area and conveyed by political authorities. Future lynx management decisions should be rooted in scientifically sound, scale-relevant information.

Gray fox (Urocyon cinereoargenteus) populations in the Midwestern USA have suffered precipitous declines in recent decades, yet they are relatively understudied. However, understanding survival and cause-specific mortality is vital for declining populations and the limited existing survival studies have been performed outside of the Midwest. We equipped 13 gray foxes in southern Illinois with GPS radio collars to investigate their survival and cause-specific mortality. We calculated the Kaplan-Meier 6- and 12-month survival rates to be 0.79 (95% CI: 0.57-1.0) and 0.53 (95% CI: 0.27-1.0), respectively. We recorded 4 mortalities: 1 disease, 1 gunshot, and 2 unknown causes. While our study has a small sample size, it contributes key information on a data-deficient mesocarnivore suffering from a population decline driven by undefined causes. We recommend further research into the survival and mortality of this elusive mesocarnivore.